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Abstract 

We present a numerical code, written in C, which can be used to simulate or to 
analyze the emission of Blazars over the entire electromagnetic spectrum. Our code 
can reproduce: Synchrotron emission, Inverse Compton emission (IC) (Thomson & 
Klein-Nishina regime). External Compton emission (EC), accretion disk variability 
using a Cellular Automata (CA) algorithm, temporal evolution of the emitting 
plasma energy distribution, flaring phenomena, light curves in the rest and in the 
observer frame (taking account for time crossing effects). 
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1 Introduction 

Blazars, are objects which emit energy in the entire electromagnetic spectrum, 
and are characterized by a variabihty which scales from long to short period. 
To understand the real nature of Blazar, we need multiwavelenght simultane- 
ous data, and we need to be able to analyze these data and to do numerical 
predictions to discriminate between the theoretical models. To the light of 
these simple but constraining remarks, we decide to write a simulator able 
to simulate the Blazar energetic and temporal behavior. A large effort has 
been produced to ensure a numerical precision enough high to use this code to 
analyze data too. This work represents just the first step of our project, which 
will have many future developments. To make our simulation as complete as 
possible, we decided to implement the simulation of disk variability too. In 
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Fig. 1. Simulation of the disk variability using CA. 

the section 2 we present the simulation of the disk, section 3 is dedicated to 
show how we can reproduce light curve both for SSC and EC model, taking 
into account for time crossing effects. 



2 Simulation of a light curve of an accretion disk, by cellular au- 
tomata 



The rapid quasi-fractal variability of Blazar is one of the main signature 
for their temporal behavior. This kind of variability is common to many 
astrophysical objects (sun spot, gal a.ctic accretion disk around bla c k hol e 



ecc....). The Cella r Automata fCA) (jMineshige. Ouchi fc Nishimori L Il994f l 



fiTakeuchi. M ineshi ge fc Negoro L Il995f ) is a fast and efficient way to simu- 
late this variability, without the need of solving differential equations. In this 
work we used CA to simulate the variability of the accretion disk emission. 
In Blazar the jet axis is close to the observer line of sight, so luminous fluxes 
are dominated by the jet emission. This led us to the understanding that the 
X-ray fluctuation, are due mostly to the emitting blob, and not to the disk. 
Therefore to have a quantitative indication about the behavior of an accretion 
disk acting around a super massive Black Hole (BH) (m > IO^Mq), we try to 
simulate the disk emission of the Seyfert galaxies, which have not signiflcative 
contribute fro m iet emission. Th e analysis of the X-ray light curves for these 



objects (Havas hida et al. 1 . 119981 ). shows a behavior similar to the 1/f noise, 
so the PSD of these light curves has almost the following form PSD oc 
with the slope a in the range 1 — 2. We than model our disk, tuning the 
input parameter, to obtain all the physical characteristic that are supposed 
to belong to an accretion disk around a super massive BH. We performed 
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Fig. 2. Simulation of EC light curves taking into account for time crossing effects 

many simulation, and in fig 1 we show just the result closest to the condition 
we believe more re alistic. To model the surface d ensity of the disk, we used 



a parametrization (ITakeuchi fc Mineshige I . Il997l ). which relates the steady 



( Shakura fc Sunvaev . Il973[ ) with the critical regime in the following way 



The parameter rj is fixed at the value of 0.5, while the value of Rcrit, has been 
varied from 10 to 90 Rs- We found that the slope of the PSD, varies with Rcru- 
The value of the slope closest to those of the Syfert galaxies is for Rcrit = 50. 
We noted how the slope change wit h the gradual mass flow rfin mA too, in good 



agreement with other simulations (IXiong.Witta fc Bao 1 . 12| 



3 Reprocessing the disk variability by the blob, a mechanism to 
explain ERC variability 



We si r aulate the temporal evolu t ion of the emitting plasma in the blob (jKataoka. J., et al. 



l2000[ ). ()Chiaberge fc Ghisellini [ Il999h . solving numerically the F okker Planc k 



equation, following the method of Chang & Cooper dch ang & C ooper , 1970l) 



and o btain light cur ves both for synchro t ron (llnoue &: Takahara 1. 1996 ). (IRvbicki fc Lightman 
979l) and IC/EC te and fc_Grindlav 1 . 1198,^^ . (|Bhimenthalfc Golddl Il970l) 



f Permer fc Schlickeiser li2002[ i emission. Than we reprocess these light curves 
taking into account for time crossing effects. Here we present a simulation in 
which the light curve of the disk (showed in fig. 1) constitutes the external 
photon field for the EC mechanism. We than assume that the Broad Line Re- 
gion (BLR) reflection of the disk luminosity (with an reflectivity x=0-0001) 
generates a quasi-isotropic photon field, and do not take into account for the 
radiation coming straightly from the disk. This light curve has a typical time 
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scale in the disk frame of about lO^s. Our simulator take into account for the 
relativistic transformation for time and frequencies, between the disk-BLR 
frame and the blob frame and between the blob and the earth. In this simu- 
lation we assume that the blob moves toward the BLR with a typical T = 10, 
so we use a beaming factor S = 10 for every the transformation. In fig 2 we 
report the result for two different values of R: lO^^cm and 5 ■ lO^^cm. The 
number of slices into which we divide the blob is such that the time crossing 
of each slice in the blob frame is of 100 s, that is much shorter than 10^/(5s 
which is the time scale variability of the disk seen by the blob frame. 



4 Conclusions 

One of the main goal of our work, was understanding the connections between 
the accretion disk variabihty, and it's reprocessing by the emitting blob. We 
observed that the luminosity fluctuations, that are present in the disk, are 
seen by the observer, if the blob radius is such that photon time crossing is 
not much larger than typical disk variability time scale. In the future we aim 
to explore the connection between the disk instability and the disk fueling 
to the light of the internal shock scenario, to simulate the synchrotron rapid 
variability too, using a self consistent physical model. 
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